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The physics of stars, their workings and their evolution, is a goldmine of problems in statistical mechanics and 
thermodynamics. We discuss many examples that illustrate the possibility of deepening student's knowledge of 
statistical mechanics by an introductory study of stars. The matter constituting the various stellar objects provides 
examples of equations of state for classical or quantal and relativistic or non-relativistic gases. Maximum entropy 

■ can be used to characterize thermodynamic and gravitational equilibrium which determines the structure of stars 
and predicts their instability above a certain mass. Contraction accompanying radiation induces either heating or 

I cooling, which explains the formation of stars above a minimum mass. The characteristics of the emitted light are 
understood from black-body radiation and more precisely from the Boltzmann-Lorentz kinetic equation for photons. 
The luminosity is governed by the transport of heat by photons from the center to the surface. Heat production by 
thermonuclear fusion is determined by microscopic balance equations. The stability of the steady state of stars is 
I— I , controlled by the interplay of thermodynamics and gravitation. 

O ' We report on a teaching experience at the EdBftTIIM^JtiMfi1E)M5|Si high-level non-specialized scientific institution 

^ in France. Its two academic years correspond to the middle years of universities in the United States; they are preceded 
Q by two years of preparation for the admittance examination, and most often followed by specialization. During the 
^Jj ] first year all students attend introductory courses in quantum mechanics and statistical physics. The course described 
^ ■ here was taken by second-year students on an optional basis. One purpose of the course was to present a systematic 
c/2 ] but elementary introduction to the physics of stars, including our Sun. For this aspect, an article by Nauenberg 
and Weisskopf^ has been a major source of inspiration. At the same time the course was used as an opportunity for 
students to strengthen their knowledge of many important topics in statistical physics, such as maximum entropy 
Q , chemical potential, equilibrium of classical and quantum gases, and microscopic transport. 

I Although the course was not intended to train future astrophysicists or professional physicists, we measured its 

(~| . success by the enthusiasm of the students and their progress in understanding physics. The appeal of astrophysics 
O is an ideal means of arousing interest in physics. Many basic laws of physics are brought into play in the workings 
^ O^ , of stars and have unexpected consequences. For instance, a striking fact is the omnipresence of quantum mechanics 
for such large objects, which is essential for explaining how a star radiates, how the nuclear plant at its center works, 
I I and why extinct stars do not collapse. A surprising feature of diffusion is also uncovered: as a result of its Brownian 
J> motion a photon emitted at the center of the Sun reaches the surface after having covered a distance of 2 x 10^° times 
the Sun's radius. 

• The necessity of relying on statistical physics is especially obvious in the study of stars, inaccessible objects known 

Qv^ ' only via the radiation we receive from them. The conditions of temperature, pressure, or density which prevail in 

■ them cannot be reproduced in our laboratories. We have therefore to use the theoretical tools of statistical mechanics 
Q^ I to deduce the properties of stars. The success of this approach is remarkable because we deal with exotic states of 
0^ ■ matter having no equivalent on Earth. 

It is also noteworthy that simple and apparently crude models are sufficient to account, often unexpectedly, for 
the main phenomena, hence offering an exceptional opportunity for getting students more familiar with qualitative 
' or semi-quantitative analyses using order of magnitude estimates. Such applications of statistical physics provide a 
i welcome balance to the more abstract approaches which emphasize its underlying mathematical rigor. Of course, 
r-| students who turn to astrophysics will have to face the full complexity of more realistic equations, which eventually 
O , lead, for instance, to a quantitative model for the Sun. At some places the course indicates how such equations can 
O ' be obtained. 

^ . Another interesting feature of this subject is that it appeals to many branches of physics, thus exhibiting the 

' unity of physics. For example, it is the only instance where the four fundamental forces of nature come into play in a 
, characteristic and spectacular manner: the stars are formed by a collapse of matter caused by gravitational attraction, 
' the light that they emit is generated by electromagnetic interactions, strong interactions provide their main source of 
energy, and weak interactions contribute in a crucial way to make their lifetime so long. The variety of the domains of 
physics involved in this course was reffected in the composition of the teaching team, including Marie-Noelle Bussac, 
Robert Mochkovitch, Michel Spiro and the authors: our background is in astrophysics and statistical mechanics as 
well as nuclear, particle and plasma physics. 

The course was delivered in 7 or 8 "blocks," each one including a formal lecture plus a class devoted to teaching 
through problems. Most of the topics covered are outlined in the following. A more detailed account of several of 
them can be found in Ref. 0. 
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This second-year course on statistical physics is not to be regarded as an advanced course, but rather as a means 
for deepening the understanding of some of the concepts aheady acquired in the first-year course. No new techniques 
are introduced. On the other hand, this course is the first contact of students with astrophysics, which dictates the 
order in which the various topics are presented. Hence, the course proceeded in a nonUnear way, which has pedagogical 
value but is not fully reflected here. Approximations often are based on qualitative arguments which partly anticipate 
the results to be established. Students may be puzzled by such a procedure before they understand that it is fruitful 
and not vitiated by circularity. After some astrophysical background (Sec. ||), the various topics are grouped into 
sections, each of which focuses on a specific aspect of statistical mechanics or thermodynamics: equations of state for 



the matter of the various types of stars (Sec. HI ) , gravitational equilibrium viewed as an application of the maximum 
entropy criterion (Se c. |rv| ), emission of light from the surface as an illustration of both black-body radiation and 
photon kinetics (Sec. M), energy transport by photons within the star (Sec. Vl), an d the use of balance equations to 



determine the production of heat in the core (Sec. VII). A synthesis is made in Sec. VIII, where further applications 
of statistical mechanics and thermodynamics to the physics of stars are suggested. Small type is used for material 
appropriate for problems. More details can be found in Refs. ^-0. 

The course begins with a brief overview off6tS^(51IM5IMMil^^^^racteristics including mass, radius, luminosity, 
and surface temperature, and how these quantities can be obtained through observations. Statistical physics enters 
immediately, for example, in the determination of the radius i? of a distant star. We first deduce the luminosity 
L, which is the total radiative power emitted by the star, from the light flux that we receive from it and from the 
estimate of its distance. We then find B from the Stefan-Boltzmann law, 

^ = '"'«''^-' "iiSsJ' <" 

where a = 5.67 x 10~^Wm~^K~^ is Stefan's constant. The surface temperature Tg of the star is obtained by 
comparing its spectrum to Planck's law (see Sec. 0A). 

A. Three types of stars 

It is useful to keep in mind various orders of magnitude corresponding to the three principal classes of stars. Our 
Sun is a typical example of the so-called main sequence. Its radius is Rq = 7 x 10^ m, its mass is M© = 2 x 10'^° kg, 
its luminosity is Lq = 3.8 x 10^^ W, and its surface temperature is Tg,© = 6000 K. Hence, its average density is 
1.4 gcm^-^, comparable to that of water on earth. The masses of all the stars lie between 0.1 Mq and 100 Mq. The 
Sun is mainly made of hydrogen, with 28% of its mass consisting of ^He nuclei and 2% of other light elements. Its 
number of protons is the order of 10^^. Heavier stars include red giants, a branch detached from the main sequence. 

A second category of stars, the white dwarfs, have masses between 0.5 and 1.4Mq, and typical radii of 5000km 
(like the Earth). Their density is of the order of lO^gcni^'^. 

The third category, the neutron stars, are even more compact objects. The density of neutron stars is the order 
of 3 X lO^'^gcm"^ and is comparable to that of the matter inside an atomic nucleus. A typical neutron star has a 
mass of 1.4 M© and a radius of 10 km. Neutron stars are observed as pulsars: we receive from them regularly spaced 
pulses of electromagnetic radiation because they rotate rapidly and radiate only in particular directions. Beyond 
3Mq, neutron stars are believed to collapse into black holes. 

B. A few queries 

One achievement of statistical physics is to explain the existence of these three very different and well-separated 
classes of stars. More precisely, we wish to understand the above orders of magnitude from simple information about 
the constituents of each type of star. 

An obvious but most remarkable property of stars is that they shine. Furthermore, in the case of the Sun, we know 
that its luminosity has remained constant over a period comparable to the age of the Earth, that is, 4.5 x 10^ years, a 
number obtained by measuring the proportions of radioactive elements in rocks. Understanding the source of energy 
for this radiation remained a challenge for 100 years. 

At this stage, simple order of magnitude estimates can be suggested to students. Assuming the Sun is merely 
undergoing some chemical combustion, how long could it keep its present luminosity? Or, assuming that its only 
source of energy is gravitational contraction, evaluate its decrease in energy since the epoch when it was a dilute gas 



(see Sec. IV). If we also assume (without justification) that this energy has been radiated at the present rate, what 
would be the age of the Sun? Such calculations were made as far back as 1854 by H. von Helmholtz, but it became 
clear when radioactive dating of rocks began 50 years later that gravitational contraction was insufhcient to explain 
why the radiation of the Sun lasts billions of years. 



2 



As soon as the huge heat associated with radioactivity was discovered, the idea emerged to look there for the 
origin of the energy radiated by the Sun. The first attempts were made in 1919 by J. Perrin and A. Eddington and a 
satisfactory explanation was given in 1937 by H. Bethe and F. von Weiszacker. However, a new question arises. The 
reactions of nuclear fusion of hydrogen into helium that take place in the central part of the Sun produce some amount 
of heat per unit time, which is exactly equal to the luminosity because the state of the Sun is stationary. However, 
such reactions are very sensitive to small changes in the temperature: they are activated by a rise, hindered by a 
decrease. Thus, if it happens at some instant that a little more power is produced in the core than what is evacuated 
by radiation from the surface, why does the internal temperature not rise, eventually resulting in an explosion of the 
Sun? Conversely if the opposite perturbation occurs, why does the Su n no t become extinct? In short, why is the 



radiation from the Sun so stable? We shall answer this question in Sec. VHl . 

We shall also wonder about the evolution of stars. How are they created? Why do they have a seemingly steady 
state in spite of their radiation? What happens after all their nuclear fuel is burnt? 

To study tai]pr'5paiMAmEBi^rS5B^S;oE*aMi«J6BSi!^,S({14MaTiC^^ of their 

evolution so as to have information about their microscopic constitution. 

A. Evolution of stars and models for their matter 

As a first approximation, the Sun is made of hydrogen. The binding energy of the hydrogen atom is 13.6 eV, while 
the molecular binding energy of H2 is lOSkcalmoP^ or cquivalently 4.5 eV. Because leV corresponds to 11600K, the 
surface temperature of 6000 K is not sufficient to dissociate the H2 molecules. However, as will be seen, the bulk of 
the Sun is much hotter, typically 2 x 10® K or 200 eV. It thus consists of completely ionized hydrogen plasma. Because 
it is globally neutral, we neglect the Coulomb interactions, and treat the Sun as a mixture of two perfect gases of 
protons and electrons with equal local densities. 

The same model holds for any star of the main sequence. It will be shown that such a star is due to the contraction 
of a hydrogen cloud due to self-gravitation. As the density and the pressure increase, the temperature and hence the 
luminosity also increase. The protostar thus formed turns into a star if the temperature of its core becomes sufficiently 
high so as to initiate nuclear reactions of fusion. The study of this evolution will rely on the model of a self-gravitating 
perfect gas of protons and electrons (Sec. ^ ). 

The emission of light in stars is maintained as long as nuclear reactions take place. However, when the nuclear fuel 
is exhausted, the production of heat in the core ceases, the internal pressure cannot be sustained, and gravitational 
contraction is resumed. For stars of masses between 0.1 and 0.5 A^©, fusion results in a plasma of He nuclei and 
electrons. As long as this gas behaves classically, it contracts rapidly, but this process stops when the electrons 
become a Fermi gas. The object becomes a white dwarf, an inert star which shines while cooling down slowly. Its 
matter can be represented by a model of a locally neutral mixture of two independent particle gases, the He nuclei 
and the electrons, with the electron gas quantum mechanical, while the He gas remains classical. This treatment is 
justified a posteriori. 

Between 0.5 and 10 M©, the stars reach a temperature which is sufficiently high so that fusion of nuclei produces 
elements heavier than ^He, such as ^^C, ^^O, ^°Ne, and ^''Mg, depending on the mass of the star. The resulting white 
dwarfs thus involve such nuclei instead of the He nuclei. The Sun will eventually become a C-0 white dwarf. 

From 10 to 100 A/q, fusion may lead to the most stable nuclei being around ^®Fe, but the temperature rises so 
much that most of the star explodes (supernova) while its core implodes, becoming a neutron star. Again neglecting 
the interactions, we can describe the constituents of such objects as a quantum neutron gas with a density comparable 
to that of nuclear matter. 

We also have to consider photon gases. Within stars photons are emitted, absorbed, and scattered by matter. 
Thermal equilibrium is thus reached locally, at a temperature imposed by the medium. 

B. Equations of state 

All the equilibrium properties of a gas of non-interacting particles are embedded in the grand potential fi(T, /i, V), 
a function of the temperature T, the chemical potential /i, and the volume V. It is given by 

^^-j dsU{e)f{e), (2) 
where M{£) is the number of single-particle states with energy less than e, and 

/(^) = ^ie-.yuT^^ (3) 
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with 1] = +1 for Bosons and 77 = —1 for Fermions. The equations of state, which relate the extensive quantities, the 
internal energy U, the particle number N , the entropy S, and the volume V , to the intensive quantities, T, fi, and the 
pressure P, are obtained from 

n^U -TS - nN ^ -PV, (4) 
dn = -SdT - Ndn - PdV. (5) 



In particular. 



N : 



-^JdeV{e), U ^ JdeV{e)ef{e) (6) 



are expressed in terms of the density of states 'D{e) = dJV/de. If a potential (here the gravitational potential) is 
included in the single-particle energies, it also contributes to /i for a given density. 

The classical limit is attained when rj can be neglected in Eq. (||), that is, when kT ^ ep, where ep is the 
function of the density defined by Af{eF) — N; for fermions it is the Fermi energy. This condition is satisfied for both 
the electrons and the protons in the Sun and for the nuclei in white dwarfs. The classical limit yields PV = NkT. 
However the nuclei do not always behave as a gas: they may crystallize in white dwarfs as studied in Problem 16.10 
of Ref. 2 by a mean-field method. The above particles are non-relativistic because their mass m satisfies mc^ 3> kT. 
We have for particles of spin s = 1/2, 

Af(.)^?l±iv/.»p.(.-|i).3^(2„.)»'^ (7) 

The electrons in white dwarfs and the neutrons of a neutron star should be regarded as cold: the opposite limit 
kT ^ £f holds. In this case {rj = —1), the Fermi factor becomes 

/(£)«0O.-e)-^(fcT)V(e-/i). (8) 
6 

We can still use the non-relativistic limit of Eq. (0) for Af{e), except for the electrons in heavy white dwarfs, where 
ep is not negligible compared to = 0.5 MeV. The above relations determine the various equations of state that 

will be needed. In particular Eqs. and (^) provide 

The uncertainty principle is reflected by the factor T? j^mV'^l^ because the momentum of each particle is order hjV^I'^ . 
The Pauli principle provides the factor iV^/'^ in ep because exclusion compels the particle number to behave as pp. 

For any non-relativistic gas (classical or quantal), the combination of Eqs. (||), (^, and (Q), that \s,M{e) cx e^l"^ , 
yields a relation between the pressure and the internal energy per unit volume: 

P=-^, (10) 

which can also be derived from kinetic theory as an exercise. 

In the ultra-relativistic limit, e ~ cp 3> mc^, which is relevant for the electron gas in heavy white dwarfs, we have 

A-W = 1^='. (") 

and the relation between pressure and internal energy becomes 

(12) 

Protostars behave as classical non-relativistic gases. However, their elementary constituents are protons and 
electrons (like in stars) only when they are sufficiently hot. Otherwise we may have H-atoms or even H2-molecules, 



4 



depending on the density and the temperature. Ionization equiHbrium associated with the reaction p + e ^ H is 
characterized by 

^J■p + fJ-c = Mh, = ^o, (13) 

where the chemical potentials are given by 

for ideal gases; C{T) is the internal partition function if the particle is composite. With Cp = Co = 2 and /cTlnCH ^ koL 
where eo is the ground-state energy of the hydrogen atom (— 13.6eV), valid for T <C |eo| /k ~ lO'^K, Eqs. (13) and 
( |l^ ) constitute the Saha formula for the ionization ratio of the H-gas. This ratio increases with the temperature, 
decreases with the density, and equals 1/2 for 

, n 6.9 X 10^ 

log ^ = 22 , (15) 

in SI units. Chemical equilibrium for 2H II2 is treated likewise, but H2 occurs only in the coldest and least dense 
clouds, at the very early stages of the formation of stars. 

For the gas of photons in thermal equilibrium which coexists with matter in all stars, we have 77 = +1 and /i = 



01 pnotons m 

in Eq. (^; Eqs. (y) and (|l2|) still hold. The density of energy, 

is generally negligible compared to that of matter, given hy U — ^NkT for kT ^ ep or by Eq. (^ for kT <ti £f- 
However at the very high temperatures attained in massive stars and in supernovae, the pressure of radiation given 
by Eqs. (|l^) and ( p^ ) becomes essential due to its rise as T^. It is responsibl e for stellar wind which expels matter. 



and for the existence of a maximum mass of about 100 for stars (see Sec. VIII) 
C. Radius of a neutron star 

It is already possible at this stage to present the basic properties of white dwarfs and neutron stars as student problems. We 
focus here on neutron stars. By estimating their density from their mass (slightly larger than Mq) and from their radius (which 
will be found to be in the 10 km range) and by using Eq. (^, we expect that ef defined by Af{£F) = N is oi the order 100 MeV. 
Even though the internal temperature is of the order of 10** K, kT/ep is as small as 10"*; the density is so large that the 
neutrons should be regarded as as if their temperature were zero. The gasjs non-relativistic because rUnC^ — 900 MeV 3> ep. 
The internal energy U is then given by Eq. and the Dressure by Eq. (y_0|); they are 4000 times larger than if the gas were 
classical. The specific heat per particle found from Eq. (p|) is 



i_dU _ TT^ k^T 



(17) 



which is weaker than that of a perfect gas by a factor of 1/3000. 

To determine the radius of the neutron star from its mass, we note that it is in gravitational equilibrium. Because the 
temperature is negligible, this condition is expressed by requiring that the total energy E = U + Eg, the sum of the internal 
and the self-gravitational energies, be a minimum as a functional of the density of particles n{r) at each point, under the 
constraint m-n / cfirn{r) = J d^r p(r) — M, the total mass of the star. From Eq. (W we obtain 

U^'-^^^[,^rHr)f^ (18) 
lOmn J 



while the gravitational energy is 



Ea^-^[d^rd^r'ppf^. (19) 
2 1 r-r' 



We use a crude model for which the density n within the star is taken as a constant.^ Then Eqs. (|l|) and yield 

o7/3 2/3^2^5/3 i 
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Eg = 



5R 



(21) 



and the minimum of E = U + Eq is attained for 



Eg = -2U, 



Gravitational equilibrium thus implies in this model a radius 



E = -U. 



/97r\2/3 



(22) 



(23) 



The more massive the star, the smaller its radius. A numerical estimate for M = Mq provides 7? — 12km, n = O.lSfm ^,ep = 
60MeV, which confirms the hypotheses on which our model is based. 

A similar theory holds for white dwarfs provided that the neutrons are replaced by electrons (for the evaluation of U) and that 
the nuclei are taken into account (for the evaluation of -Eg)- Their radius and density are found as function of their mass, with 
numerical values in agreement with the data of Sec. IIA. 



D. Stability problems for neutron stars and white dwarfs 

For more massive neutron stars than considered above, the Fermi energy becomes significant compared to rrinC^ and the 
relativistic form e = (m^c'* + c^p^)^^^ of the kinetic energy should be taken into account. This form produces a crossover from 
Eq. (0) to Eq. (y) for J\f{e), and from Eqs. (^ or (§o|) to 

rr ^3 r , ^l4/3 3^^^ hcM*^' 1 

U = -TTT d r \n[r)\ ' ~ tttt-;:; 24 

8^1/3 J L ^ 2"/3 7r2/3 ^4/3 i? ^ ^ 

for U. Hence, for sufficiently small R, E = U + Eg is linear in , with a coefficient which becomes negative for 

Thus, neutron stars heavier than this limit would collapse because the Fermi pressure is insufficient to resist gravitational 
attraction. 

As shown later, this calculation is irrelevant because another effect prevents the existence of such massive neutron stars. 
However, for white dwarfs, the corresponding calculation yields the Chandrasekhar limit mass beyond which white dwarfs cannot 
exist. The occurrence of this maximum mass implies an upper bound on the luminosity of white dwarfs. The Chandrasekhar 
mass is found to be 1.7 Mq by a crude estimate analogous to Eq. (|2^), and to be 1.4 Mq by a more precise calculation taking 
into account the complete form of the relativistic energy e and the variation of the density. 

Returning to neutron stars, we note that Eq. (^) implies that ep = ^GMnin/ R- Hence, the total energy of a neutron at the 
surface of the star, e — GMrrin/ R, is negative because e < ef- No neutron can escape from the star because the escape velocity. 



1/2 / \ 2/3 



( 2Gm \ 

I Mq 



(26) 



is larger than the Fermi velocity. The velocity in Eq. (|2^) reaches c for M = 3Mq , suggesting that photons themselves cannot 
escape if M > 3Mq. Thus there exists a maximum mass for neutron stars, beyond which such objects become black holes 
where light is trapped by gravity. This mass is smaller than the mass of collapse in Eq. (^5|) which is therefore irrelevant for 
neutron stars (but not for white dwarfs). The present argument is faulty because it uses the non-relativistic approximation for 
the energy of neutrons although £f equals jm-aC? when the velocity in Eq. (^^ reaches c; furthermore general relativity should 
be used to describe such strong gravitational fields. However, the conclusion turns out to be correct.^ 
A neutron decays as 

n ^ p + e + 77, (27) 

with an average lifetime of 15 minutes. How can a neutron star be stable against such a /9-decay? Actually the electrons and 
protons created by the reaction (|27| ) can in turn react as 

p + e ^ n + i/, (28) 
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so that a stationary regime sets up where the star contains, beside A^n neutrons, a number A'p — of protons and electrons. 
The neutrinos practicaUy do not interact with matter and are radiated, only carrying out a small amount of energy. If we 
forget about them, we determine A^c by the chemical equilibrium condition 

/in = Mp + Mc- (29) 

The chemical potentials include contributions mc^ which can be omitted because (mn — mp)c^ = 1.3 MeV and rricC^ = 0.5 MeV 
are negligible compared to the Fermi energy ep of the neutrons. They also involve contributions from gravity which are almost 
the same on both sides of Eq. (|2E|). Assuming that the protons and the electrons constitute "cold" Fermi gases, and assuming 
the electrons to be ultra- relativistic — hypotheses to be checked in the end — we note that A^p = A^e implies that the Fermi 
momenta Pp = pp are equal and hence that Ep <^ £p. Thus A'e results from ep = Ep, which together with Eqs. (H) and ( |l l| ) 
yields 

^-^(^X^ (30) 

which is 2 X 10~^ for M = Mq. The presence of this small number of electrons and protons is sufficient to prevent the /3-decay 
of the neutrons. 



IV. SELF-GRAVITATING OBJECTS: MAXIMIZING THE ENTROPY 
A. The relation between mass, radius and temperature 

Like any thermodynamic equilibrium, gravitational equilibrium is expressed by looking for the maximum of the 
entropy, subject to constraints on conserved quantities. The existence of radiation within a star causes a negligible 
departure from equilibrium (except for very massive objects), so that we have to write that S is maximum for a fixed 



value of the total energy E = U + Eq- We begin with the same crude approximation^ as in Sec. [II C, neglecting the 
variations of density and temperature within the star or protostar. We thus have to maximize 

S-i3{U + Eg) (31) 

with respect to U and to the radius, where p is the lagrangian multiplier associated with the constraint on E. The 
first condition implies that = [3. If we use Eq. (pT| ) and vary the above with respect to R, we obtain 

where P ~ ^dU/dV is the average thermodynamic pressure. For a given mass M and a given temperature T, the 
radius R is determined by Eq. (p2|), where P is a function of T and the density. 

The long-range nature of the gravitational forces implies a non-extensivity of the matter of stars. Doubling the 
mass of a star for a given temperature does not double its volume. In particular, the volume of a neutron star or a 
white dwarf is inversely proportional to its mass (see Eq. (p3|)). 



From the microscopic viewpoint, the gravitational attraction at short interparticle separations produces another pathology. The 
Boltzmann-Gibbs distribution e^^^'''^ contains in classical statistical mechanics a factor exp ^ Grn^ /rijkTl for particles with 

U>j -I 

masses m and relative distances rtj. This factor becomes so large at short distances that the partition function diverges, and 
hence, there exists no equilibrium in classical statistical mechanics. The expressions or (Ell) for the gravitational enere'v 
implicitly rely on a mean-field approximation and thus disregard the correlations wlucn invallaate the classical calculatii 



tion. 



However, as seen in Sees. |l|c and HID, the Pauli principle cures this difficulty by providing a short-range cut-off ft/pr, and 
thus prevents the occurrence of the collapse, at least for not too massive stars. 



The condition ( p2[ ) characterizes a stable equilibrium provided 5 is a maximum, that is, d S 
d [T-'^{dU + PdV)] < for d£' = and dS = 0. This condition is expressed by 



dP_ 
dV 



d^EG _AEg _^AP 
s dV^ ~ 9V^~ 3V 



Hence, the gravitational equilibrium is stable only if the adiabatic compressibility ks = ~ {P/V){dV/dP)\g is smaller 
than 3/4. This condition is satisfied for non-relativistic gases, whether classical or quantal (for which ks — 3/5), 
but not for ultra- relativistic Fermi gases (for which /tg = 3/4). Such a relativistic instability is responsible for the 
gravitational collapse of massive objects, in particular for the implosion of a supernova leading to a neutron star. 
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For non-relativistic gases, Eqs. (10) and (p2f) imply that 



E^U + Eg^^Eg^-U ^-^PV, (34) 

whereas for ultra-relativistic gases, Eq. ( p^ yields the total energy 

E = 3PV + Eg^O, (35) 

which thus vanishes when the star becomes unstable. 

If the stellar object radiates, and if the corresponding loss of energy given by Eq. (|^) is not compensated for by a 
production of heat in the core through nuclear fusion, the decrease \dE\ in energy produces a shift of the gravitational 
equilibrium. For non-relativistic gases, the resulting changes in the various thermodynamic quantities are obtained 
from Eqs. (^)-(^, by noting that the Massieu potential per particle 4'(/9, V) — {k^^S — l3U)/N is a function of the 
dimensionless variable /3~il/2/3^-2/3 ;j-2^^ (This result can alternatively be found from Eq. (^J).) Hence, using the 
fact that the first and second partial derivatives of ^' with respect to /? yield U and the constant volume specific heat 
Cv — dU/dT\v, one can find as an exercise: 

dE dR IdP dU TdS CvdT 

(36) 



U R A P U U 2CvT~U' 

The temperature may increase or decrease, depending on the sign of 2CyT — U . The entropy always decreases^ even 
when the temperature rises, due to the contraction which accompanies the loss of energy. A larger increase of entropy 
takes place in the outside world, because an amount \dE\ of energy is radiated toward regions with temperature lower 
than T. 

We recover from the above results those of Sees. QIC and |i|d on neutron stars and non-relativistic white dwarfs. 



In these cases, CyT is neghgible compared to U , and radiation induces coohng. 

B. Evolution of a hydrogen cloud 

A hydrogen cloud behaves as a perfect gas, so that 

-E^U ^ ^NkT = CvT, (37) 

where TV is the total particle number. Hence, a loss of energy \dE\ is accompanied by heating. Gravitational forces 
thus induce an effective "specific heat" dE/dT = —^Nk which is negative. In contrast, for extensive states of matter 
without gravity, the concavity of entropy implies that ^' is convex and hence that the true specific heat dU / dT is 
positive. 

As the cloud radiates and shrinks, the temperature rises, and the heating becomes more and more efficient, 
especially if the cloud has a large mass. Once the ionization temperature given by Eq. ( p^ ) is reached, the particle 
number is = 2M / rrip and the condition (^) expresses the radius as 

R=^^. (38) 
lOk T ^ ^ 

This relation holds for protostars as well as for stars of the main sequence for which the electron gas is classical, with 

We denote hy n — M/mpV the common density of electrons and protons (see Figs. 1 and 2). For the Sun, Eq. (|3^ ) 
yields an average internal temperature of 2 x 10^ K. 

A star of the main sequence is in a stationary regime, where the luminosity L is exactly compensated for by the 
power Q create d in the core by nuclear reactions. Its temperature remains fixed at a value about 10'' K or 1 KcV in 



the core (Sees. VH and VIII ) . Its radius has the value corresponding to Eq. ( |38| ) and the energy E (excluding the 



m 



V 



nuclear mass energy) remains fixed at the value —3MkT/ 

For a protostar the radius decreases with the total energy as shown by Eq. ( |38[ ) , while the temperature rises until 
it reaches (for M > 0.08 M©) the ignition value for H, where the contraction stops. Clouds with masses M < 0.08 
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do not result in stars; as seen below, their temperature never attains the ignition value. The cloud results in an 
aborted star, a brown dwarf, or for smaller masses an object hke Jupiter (M = 10"^ Mq). 

Once a sizeable proportion of the protons in the core has been transformed into helium, its fusion can no longer 
take place. Radiation is no longer compensated for by heat production, which entails anew a decrease of hence 
shrinking and heating. If the star is sufficiently massive [M > 0.5 Mq), the temperature may then reach the higher 
ignition value for fusion of He into C and O. For still more massive and hotter stars (M > 8Mq), fusion can lead to 
heavier elements. 

During these processes, the left-hand side of Eq. (|39|), proportional to T'^/^M"^ because n oc MR~^ and R cx 
MT^^, increases until the condition ( |39| ) is violated. The electrons gradually enter the quantum regime, described 
by Eqs. while the classical approximation still holds for the residual hydrogen and for the nuclei produced. 

According to Eq. (^), while the radius always decreases when the either aborted or extinct star radiates its energy 
without thermonuclear production, the temperature first rises, reaches a maximum, and then decreases: while dT ~ 
\dE\ jCy in the classical regime of Eq. (p9|), we find from Eq. ( p6| ) that dT = — \dE\ jCy in the Fermi gas limit, with 
Cv oc T. 

More precisely, the relations between the density, the temperature, and the energy for a given mass can be found by using 
Eqs. (g), (^, and and Eq — —2U. They result from the elimination of the electronic chemical potential fic = kTa between 

_ iV, (2mefcr)3/^ 

= IT = 12^2 fi3 h/2{a), (40) 

and 

E _ GM^ _ iVe /5/2(q) , iVnuc 



3kT WRkT 5 13/2 (q) 



(41) 



The quantities A^c and A^nuc are the total numbers of electrons and nuclei, respectively, and 

/o cosh^\{x-a) ^ ^Jo e^-'" + 1 



The maximum temperature is obtained by cancelling 

2CvT -U = 3kT 



4/5/2(0) 73/2(0) ^^^ , A^n 
5/3/2(0) h/2{a) 



(43) 



which provides a ~ 4.44; then by using Eqs. (0) and (^), we find 



r^ax = 2.5xl0'(-^)*/^K. (44) 

lvl( 





The numerical values were estimated for an aborted star where A^nuc — Nc — M/nip, but this behavior is general. As expected, 
higher maximum temperatures are reached for more massive stars. When the star cools sufficiently, it becomes a white dwarf, 
described by letting o ^ 00 in Eqs. (^0|)-(^2|). In particular E and R end up at a minimum value given by Eqs. (po|)-(P3[), 

where rUn^^ should be replaced by mclmpA/ZY^^'^ and Z and A are the average charge and mass numbers of the nuclei. It is 
instructive to represent the above results as graphs in the density-temperature plane (Fig. 1) and in the density-energy plane 
(Fig. 2). 

A numerical estimate of the thermal plus gravitational energy of the Sun from Eqs. (^^ and (jsij) shows that it is 1000 times 
smaller than the radiative energy released during the whole life of the Sun (for 10^" years at the rate of Lq = 3.8 x 10^® W) . 
Thus the energy of gravitational contraction has been negligible compared to the nuclear energy released and radiated; it has 
only contributed to bringing the core of the Sun to the ignition temperature for the fusion of hydrogen. However, once the 
nuclear reactions stop, a new contraction will produce an overheating and a huge increase of luminosity. The amount of energy 
released in this process can be estimated from the final energy of the Sun as a white dwarf, given by Eq. ( jZlj ) for a — > 00. It 
is found to be order of one tenth of the total radiated energy. Thus, in the stage following exhaustion of the nuclear fuel, the 
radiative energy feeded by gravitation only will be a significant fraction of the thermonuclear energy. 



C. Local density 

A more precise theory requires us to account for the non-uniformity of the mass density p{r) and the temperature 
r(r) within the stellar object. We assume that the total angular momentum vanishes, so that gravitational equilibrium 
ensures spherical symmetry. The entropy S and the internal energy U are now sums of contributions Si and Ui from 
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small volume elements uji, and Eq is given by Eq. (pj|). Equilibrium is expressed by maximizing Eq. (|3^) with respect 
to variations Sui and to deformations which displace matter from r to r + Sr, where dr is a function of r. The volumes 
uji are changed according to Suji = UiV ■ (5r, and the variation of Eq is 



5Eg = j (Fv p{r) VVF • ^r, (45) 

where W{r) is the gravitational potential 



r 



r'2 dr 







/•OO 

'p{r') - AttG / r'dr'p{r'). 

J r 



The variation of S is 



- E [ t + t ^-«] - E t V . ^r. (47) 

i i 

Thermodynamic equilibrium thus implies for independent variations 6ui that Tir) is uniform and equal to I3~^ . It 
also implies for arbitrary deformations that 

VP = pVW (48) 

at each point. We thus recover the equation of hydrostatics as a consequence of the laws of thermodynamics. 

Because P is a function of T and the density, Eq. ( ^ ) is an integro-differential equation for the mass density, to 
be solved with p{oo) = and / d^r p{r) = M. It is equivalent to the differential equation 

VP 

V = -AnGp. (49) 

P 

A consequence of Eq. ( |4^ ) is the virial theorem 

3ld^rP=-EG. (50) 



It can be proved directly by writing that the thermodynamic potential ( pT| ) is stationary under a dilation = re, 



Sui — —PiSuji, for which 5U = — 3/d^rPe, and 6Eq, — —Eqe. The approach of Sees. IVA and IVB was based on 
the virial theorem in the crude form of Eq. (p2[). 



Actually, stars are not in thermal equilibrium because they radiate and because heat is produced in their core. The temperature 
T(r) decreases from the center to the surface. The time-scales involved for thermal exchanges are large (it takes 10^ years for a 
photon to propagate through the Sun), whereas gravitational equilibrium is reached very rapidly due to the long-range of the 
forces (it takes only a few hours for sound waves to propagate across the Sun). Likewise, electric neutrality is locally ensured 
by the long-range Coulomb forces. In the above calculation, the variations Sui and Scui thus correspond to mechanical energy 
only, without heat transfer nor diffusion. They are thus not independent, but satisfy 5ui = —PiSui. We thereby obtain no 
condition on T{r), but Eq. (^) stiU holds. 

In stars which have already transformed a significant part of their hydrogen into helium and heavier elements, equilibrium is 
also not reached regarding the variation in space of the proportion of elements. This situation occurs in particular for red 
giants and for supernovae. In the gravitational potential W , the density of nuclei with masses would vary in equilibrium 
as &qi[—miW (r) / kT], whereas nuclei remain in stars at the place where they were formed because their diffusion is a slow 
process. (Convection can however raise heavier elements up to the surface.) The above reasoning, in which matter is deformed 
as a whole, applies here and leads again to Eq. (^8|). 

As an example, for a not too massive white dwarf (with M < O.SM©), the thermodynamics is governed by the electron gas 
which is non-relativistic. The nuclei, which for such stars are H or He, only contribute to the gravitation. Equations ^ and 
(^) hold for the local density of electrons ne{r), the internal energy per electron, and the pressure P{r), which are related to 
one another by the local Fermi energy eF(r). The mass density equals 

P(0 = ^n.(.) = ^^|5g^[.p(r)r, (51) 
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where A and Z denotes the average mass and charge numbers of the nuclei, because on average one electron drags 1/Z nuclei. 

-PdV, that |ncdP is the differential of the energy per particle given by 



On the other hand, we find from U = ^PV and dU ■ 



Eq. 



so that 



ivP=^V£F. 

p nipA 



(52) 



The differential equation 



for p{r) can thus be rewritten as 



r 



MO] 



3/2 



(53) 



in terms of the dimensionless variables 



^(0 = 



£F(r) 
eF(0) 



p(r-) 



P(0) 



2/3 



£f(0) 

p(o) ■ 



(54) 



The solution of Eq. (p^, the Lane-Emden equation of index 3/2, with boundary conditions f{0) = 1 and ip'iO) — 0, is a function 
which decreases from 1 to 0, a value attained for = 3.65. At this value <p{(,) has a finite derivative, £.iip'(£,i) = —2.71, and 
beyond it vanishes. The parameter a, and consequently £f(0) and p(0), are determined by expressing the total mass as 



M = / d^rp{r) = 47raV(0) j ^'c^d^ MOf^^ 
t2 N 3 



9n^ 



2mcGa 



CiV'(6)| 



(55) 



where we have used successively Eqs. (^ij), (^l|), and (^3|). 

We see that a universal function </p(5) describes the density of all white dwarfs with masses M < 0.5 Mq. A sharp edge occurs 
at the radius _Ri — which scales with the mass as M^^^ , and which equals 1.2 times the radius provided by the crude 
approximation of Eq. (p3|). Near _Ri, the density decreases as {R\ — r)^''^. The central density p(0) also scales as R^^^ . Its 
value is larger than the one found in Sec. |ll|c by a factor of 3.3. For white dwarfs with masses M > Mq, the large value of 
p(0) implies the failure of the non-relativistic approximation, and e — /2m should be replaced by {m?c^ + (P'p^Y^^ for the 
electrons in the central part of the star. When M approaches the Chandrasekhar mass, the ultra-relativistic approximation of 
Eq. ( |ll[ ) may be used in that region. 

Because for such stars, P depends only on the density through Eq. (p2|), the hydrostatic equation (Uq) can be integrated as 



£F(r) + ^Vy(r) = p. 



(56) 



The integration constant p is interpreted as the total electronic chemical potential, including the contribution of the effective 
gravitational potential seen by an electron together with the nuclei that it drags. Indeed, equilibrium implies that the chemical 
potential is uniform, at least when the temperature is also uniform or as here negligible; otherwise, its gradient is obtained 
from Eq. (fig) and from the Gibbs-Duhem relation as — sVT, where s is the entropy per particle. 



V. THE LIGHT EMITTED BY A STAR: PHOTON KINETICS AND BLACK-BODY RADIATION 
A. A star as a black body 

The radiation of a body is characterized by its luminance Li,{6), defined as follows. The power 5w radiated by a 
surface element 5S in a solid angle 5liJ around a direction making an angle 9 with the normal to the surface and in a 
range of frequency v,v + 5v equals 

6w = L^{e)cose5S5uj5v. (57) 

The main features of light emission by a star can be understood by regarding its surface as a black body, that is, as a 
material which is a perfect absorber. Standard thermodynamic reasoning shows that the luminance Ll{T) of a black 
body at temperature T is related to the spectral energy density of a photon gas in equilibrium in an enclosure. From 
Eqs. (^, and ( pT| ) we can obtain Planck's law 
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Photometry shows that the spectrum of the light that we receive from stars can be fitted to Eq. (58), the 
parameter T being identified with the surface temperature Tg of the star (5800 K for the Sun). As a function of the 
frequency, has a maximum at /ii^max — 2.82 kT (Wien's law), so that the color of a star gives us direct information 
about its surface temperature. One classification of stars is actually based upon their color. 

The total luminosity L of a star is obtained by integrating Eqs. ( |57| ) and ([s^), which yields Eq. (|l|). On the 
Hertzsprung-Russell diagram,** stars are plotted according to their color and luminosity, which exhibits regularities 
that will be explained in Sec. VID. 

The fact that does not depend on is Lambert's law. It expresses the fact that the power 5w received by an 
observer is proportional, not to the surface 5S, but to the apparent surface 6S cos9 of the emitter. The angle 6 is not 
accessible for distant stars, which are always seen as points. For the Sun, 9 increases from at the center of the solar 
disk to 7r/2 at its edge. We shall discuss the validity of Lambert's law in Sec. ^C. 

B. Kinetics of photons in matter 

A microscopic theory is needed to explain the above results and to obtain corrections. We must take into account 
the structure of the matter of the star, and in particular, the fact that its surface is ill-defined. The light comes from 
a superficial shell which is partly transparent, the photosphere. In this region the temperature T(r)and the density 
p{r) decrease with r and we need to understand the meaning of the quantity Tg which enters Eqs. (|5^) and (|l|). The 
propagation and thermalization of particles moving in a gaseous medium, such as photons in the photosphere, is 
described by the Boltzmann-Lorcntz kinetic equation for the particle density f{r,p,t). This equation has the form 

^+v.V/=X(/). (59) 



For photons with energy e — hv = cp, it is convenient to rewrite Eq. ( |59| ) in terms of 

A,(r,u,t) = ^/(r,p,t), (60) 

where u = p/p is a unit vector in the direction of propagation. When r lies just outside the star, A^, can be identified 
with the luminance entering Eq. (|57|). The kinetic equation then reads 

+ CU • VAi. = Jc + ^'sp + Jst - ^a, (61) 

where the right-hand side depends on the distribution function Ai, and on the local properties of matter. The term 
is the collision term describing elastic or inelastic scattering of photons, Xsp and 2st are the source terms associated 
with spontaneous and stimulated emission, and is the absorption term. At each point, matter is in local thermal 
equilibrium at temperature T{r), but the photons are far from equilibrium, especially near the surface where their 
distribution is extremely anisotropic as they are escaping. In the photosphere, T lies below the ionization temperature 
and the matter mainly contains atoms which govern the dominant processes, absorption and emission. Inside the star 
where matter is completely ionized, the dominant process is the scattering by electrons. 

We focus here on absorption and emission. The corresponding three terms in Eq. (^l]) are associated with the 
reactions 7 -I- Ai ^ A2, where the states Ai and A2 of an atom have energies ei and £2 = ei + hv. Because the 
photosphere is in thermal equilibrium, the densities ni and n2 of the two species satisfy 712 = nie'^^/'^'^ . The average 
number {N^) of photons in a single mode with frequency v is 

(iV.) = \h^f = ^A.. (62) 

Denoting by (7^. the cross-section for absorption of a single photon by a single atom, the flux of photons reaching 
this atom is c/ and the absorption term in Eq. (|6^) is therefore X^ = cCaniA^. Probabilities of opposite elementary 
processes are equal. Hence, the simulated emission term is Xjt = c(Tan2A,y, and the spontaneous emission term Xjp is 
obtained for each photon mode by dividing the contribution of stimulated emission by the number of photons. (A 
similar relation between two inverse processes holds for inelastic scattering.) With these considerations and Eq. (|6^), 
Eq. ( |6l| ) becomes 

-^+U- VA„ = Vda 712 A, + ^ -71iA, . 63 

c at ^ — ' 
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The summation is carried over the various processes which can occur due to the existence of various constituents in 
the photosphere. Using the relation ni = 712 e^'^^'^'^ , we can rewrite Eq. ( |63|) as 

i^+u-VA, -fc,(ilJ-A,), (64) 
c ot 

where the coefficient k^, is defined by 

^ CTa(rii - 7^2). (65) 
The quantity can be interpreted as a mean free path for the photons with frequency ly. The quantity on the 



right-hand side of Eq. (|6j) is the same as in Eq. (58), evaluated at the local temperature T(r) of the medium. 

As expected from equilibrium statistical mechanics, the right-hand side of the kinetic equation vanishes when 
the photons are thermalized with matter, their distribution then being the same as within an enclosure at temper- 
ature r(r). The temperature T(r) appears not only through L^, but also through the factor fcj^, which moreover is 
proportional to the local density p{r) and depends on the frequency through the cross-sections a^.. 

Regarding the surface of a star as a plane, the quantities and are functions of the altitude z. In a stationary 
regime, Ai,, a function of z and of the angle 9 between the direction of propagation and the 2;-axis, should be obtained 
by solving 

^ = ^ALl~K). (66) 
dz cos 

If the surface of the star (at z = 0) were sharp, and if the temperature and the density were uniform for 2: < 0, the 
solution of Eq. (^) would be L° = everywhere for cos^ > 0, and would satisfy A^ = for 2 > and cos 9 < (no 
incoming light from outside the star), and go to for z — > —00 and cos 9 < (light propagating inwards inside the 
star). The star would thus radiate like a black body as indicated in Sec. ^A. 

To solve Eq. ( |66| ) for the realistic case of an inhomogeneous photosphere, we define the optical depth as 

/>oo 

C=/ dz'k,{z'). (67) 



This quantity is dimensionless and depends on the frequency and on the geometrical depth z. It vanishes outside 
the star, increases when the altitude z decreases, and tends to infinity deep within the photosphere (2 — > —00). The 
optical depth characterizes the penetration of light in the photosphere, which is governed by the cumulative effect on 
photons of matter lying above 2; an incoming beam with cos 6' < would be attenuated as e~^/l'^°''^l as it propagates 
inward. 

Inverting Eq. ( [67[ ) gives 2 as function of ( for a given frequency. Thus the luminance can be regarded as a 
function of and the kinetic equation (|66|) for A^{(,9) reads 



dC, cos 9 



(68) 



Equation ( |6q ) should be supplemented with boundary conditions expressing for cos 6* > 0, the boundedness of K^, deep 
within the photosphere (C — > 00), and for cos6' < 0, the absence of incoming light outside the star (A^, — for ^ = 0). 
The solution is 

KiCe)^ r ^c-^'/^-^'^LliC + C), (O<0<7r/2), (69) 







cos 9 



and 



Jo \cos9\ 



(70) 



Using Eq. (p3) outside the star where C = and where A^ — L^, we find the luminance of the star 



LA0)=l -^e-^/^°^'L'ATiz)], (71) 
cost^ 
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where L° should be evaluated from Eq. ( |58| ) at an altitude z which is related to the optical depth C and to the 
frequency by Eqs. ( |65|) and (|6^). 



C. Departures from black-body radiation 

The combined effect of emission and absorption has resulted in the expression (71) for the luminance of the star, 
which for fixed Q and v appears as a weighted superposition of black-body radiation associated with successive layers. 
Contributions from deeper and deeper layers involve larger and larger values of because the temperature increases 
with the depth, but they are weighted by e~'»/'^°''^, a decreasing function of the optical depth. The region which 
contributes to the radiation is therefore the one in which C rises from to some number of order one. 

The photosphere thus defined is thin. For the Sun, the optical depth C rises from to 1 in a region of only 
20 km thickness. It reaches 10 when the altitude decreases by 500 km, small in comparison to Rq — 700 000 km. 
Because the photons that we receive all come from nearly the same altitude, the edge of the Sun appears optically 
very sharp, although the density remains significant in the chromosphere, a shell with thickness 1500 km which lies 
above the photosphere. Moreover, because the photosphere is thin, its temperature is nearly uniform, which explains 
why the black-body approximation of Sec. is fairly good, both for the total luminosity and for the spectrum of 
the emitted light. Indeed the temperature decreases in the pho tosphere of the Sun by 7% for an increase by 50 km of 
the altitude, and the parameter Tg, which in Eqs. (|l|) and ( p8|) is fitted to the observed radiation, corresponds to an 
average temperature in the photosphere. 

A better approximation consists of linearizing T(z) and L^{T) in the photosphere. Moreover, we neglect the 
variation in the photosphere of k^{z) defined by Eq. (Ba). Equation ([7l|) then provides, using C — ~ k^^z, 



( 

cos 6 



dz 



(72) 



For given and 9, the luminance L^{6) of a star is therefore the same as that of a black body at the temperature 
which prevails at the altitude z = — cos6/k^ such that the optical depth ( equals cos6'. 

The dependence of L^{0) on the frequency arises both directly in and through the temperature T{z = 
— cos0/k^). The contribution to k^ of inelastic processes is fairly constant. However, we see from Eq. ( |65| ) that 
atomic transitions correspond to sharp resonance peaks in the absorption cross-sections and hence in k,^, with widths 
inversely proportional to the lifetime of the excited level. Around such a peak, the depth \z\ — cos9/k^ has a narrow 
minimum, and Ll{T) should therefore be evaluated at an altitude where the temperature is significantly lower than 
outside the resonance. Thus, for each atomic resonance, the luminance has a sharp minimum as function of the 
frequency, and the spectrum displays dark lines associated with an increased absorption. The intensity of light at the 
resonance frequency is determined by the temperature of the layer involved. 

A converse effect also exists. In the chromosphere and in the corona which extends beyond the altitude of 
2000km, the temperature strongly rises reaching 10^ K. The density is very low, which strongly reduces k^. However 
at a resonance, fc^ has a significant value and the optical depth is no longer negligible. The very high temperatures 
involved in L1{T) thus produces bright lines associated with an increased emission, compared to the average Planck 
spectrum associated with the temperature of 5800 K. 

The dependence of L^{9) on the angle 9 arises from the factor cos6' entering the altitude z = — cos 9/ ki, in 
Eq. (|7^). The light that we receive from the center of the solar disk corresponds to [T{z — — l/k,^)], whereas the 
edges send us light associated with LjJ [T'lO)]. Thus the periphery of the solar disk appears both less bright and 
more red than its center, because the effective surface temperature which characterizes the radiation is lower there. 
Photographs show these features. If the Sun behaved as a black body with uniform temperature, Lambert's law would 
be satisfied, and the Sun would appear to us as a uniformly bright, seemingly flat disk. 

When the radiation of the photosphere is occulted by a total eclipse, the corona becomes visible. Its luminance is weak in spite 
of its high temperature, contrary to what happens for a black body. This can be understood as a consequence of KirchhofF's 
law, which expresses the luminance of a body as the product of the luminance L° (T) of a black body at the same temperature 
by the absorption coefficient. This coefficient is here very small because the low density of the corona makes it transparent. A 
quantitative approach relies on the kinetic equations (pl|). Disregarding the scattering of photons but now taking into account 
the spherical geometry of the Sun, we can find as an exercise the full stationary solution of Eq. (p4|). Hence we derive the 
luminance (x) which can be observed along a line passing at the distance x from the centre of the Sun. For x < Rq we 
recover the result of Eqs. ( [n[ ) or (^^ with x = _RQsin^. For the outer atmosphere x > Rq which is transparent, we have 
k^RQ <^ 1. In this limit the ruminance of the corona is found to be 

U{x)^2 Z^ fc. (r)L° [r(r)], (73) 
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a quantity much smaller than the black-body radiation at the temperature of the corona, but still following Planck's spectral 
distribution at this temperature. 



Quantitative studies of L^{9) as function of v and provide extensive information on the elements constituting 
the external shells and on the variations of the temperature and the density with the altitude. 

TheYJfetB-Affli^^a'lB;^^ mA(M§EiifissWSSHK>f lfeme&sPfi¥dJS£a^^h4^^^fitSRaW^ii^d[AmMa1^^^ 

in the c ore e ither by gravitational contraction (protostars) , or by thermonuclear reactions (stars in the main sequence, 



see Sec. VII), or it is stored there due to earlier evolution (white dwarfs). The heat flow is mainly induced by thermal 
conduction in white dwarfs and by convection in red giants. Convection is also important in the upper shell of the 
Sun (1/4 of its radius). However, the dominant process of energy transport within protostars and stars of the main 
sequence is radiative transfer by photons which are emitted, absorbed, and scattered by matter. 

A. Transport equations 

At the microscopic level, the dynamics of photons inside a star is governed by the kinetic equation (^9|) or ( |6l| ) 
as in the photosphere. However, not only the matter but also the photon gas are now in local equilibrium due to 
the effects of emission, absorption and inelastic scattering processes. We therefore can use the macroscopic equations 
of non-equilibrium thermodynamics. The distribution Ai,(r, u) of photons defined by Eq. ( |60| ) is close to the black- 
body distribution [T{r)] of Eq. (|58|), but the current density of energy J(r) for photons, which is the dominant 
contribution to the heat flux in the considered regime, originates from the deviation — from equilibrium. It is 
given at each time by 

J(r) = Jd^pcuf{r,p)cp^ j d?udvuky{v,u). (74) 

Conservation of energy is expressed by the balance equation 

3u 

^ + V-J-9, (75) 



where u is the total energy density, mainly carried by matter and evaluated in Sec. HIB. The source term q{r) on the 
right-hand side describes the heat production by thermonuclear reactions. 

Thermal conduction is characterized by the heat conductivity, the ratio between the heat current J and — VT. 
For photons, it is customary to replace VT by the gradient of the energy density of photons, 

u^{T)= Jd^pf{r,p)cp~^Jd^udiyLl{r) = ^aT^. (76) 

This replacement defines a transport coefficient, the opacity k, through 

3 = -^Vu^. (77) 
3pK 

The mass density p of the medium has been introduced in Eq. (|77|), because for a given constitution of matter, the 
fiux of photons generated by a fixed thermal gradient is inversely proportional to the density of particles on which 
they scatter. The opacity is measured in m^ kg^^, and the quantity 

£ = — (78) 

pK 

can be interpreted as the mean free path of photons between two successive collisions. 

Equations ([TS]), (|76|), and (^^ are sufficient to determine the transport of energy in terms of the opacity k. 
The value of k can be evaluated by solving the Boltzmann-Lorentz equation (|6l|). Within a star, when ionization 
is complete, the main radiative process is the elastic Thomson scattering of photons by electrons, for which the 
cross-section, 

duj Att 6 ATreomcC^ 

is isotropic. From the definition of a cross-section, the elastic collision term in the kinetic equation (|6l|) is obtained as 
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Ir. = cn, I d?W ^ [A,(r, u') - A,(r, u)] 



(80) 



This collision term dominates the right-hand side of Eq. ( |6l|) . However the remaining absorption and emission terms 
given by E g. (34| ) are not irrelevant: they ensure thermalization, and allow us to replace by in the left hand 
side of Eq. (^). Multiplication by u and integration over the direction of u of Eq. (pH ) with the right hand side (|80| ) 
then yields, in a stationary regime, 



j d^uu (u • Vi") -JicCTh j d^uu A,,. 
Using Eqs. ( [7^ and (|7^) and integrating Eq. (|T]) over v, we find 



which proves microscopically the empirical equation ( [77[ ) and provides 

«=^^4xl0-2m2kg-\ 



(81) 



(82) 



(83) 



B. The random walk of photons 

Let us first disregard the variations of density and temperature in the medium. From Eqs. ( [75|) and ( [ttI ) wc find 
for u~f, the density of energy of photons, the heat diffusion equation 



du^ c 
dt 3pK 



V u-y — q. 



(84) 



According to Eq. (p4|), a unit pulse of energy which is given to the photons at the point r = at time < = diffuses 
for i > as 



3pKr 
Act 



(85) 



The perturbation spreads as (r^) — {2ct/pK)t. 

Taking for p the average density of the Sun, 1.4 x 10'^ kgm~'^, we find that it takes 1500 years for (r^) to reach Rq. 
If the matter were transparent, the photons created at the center would reach the surface after only 2 seconds. Their 
scattering by the electrons slows them down enormously. The opacity for the photons of stellar matter is therefore 
essential for preventing the energy confined in the core from escaping immediately. In contrast, matter is transparent 
for neutrinos, so that most of the available energy in supernovae is suddenly evacuated in the form of neutrinos. 

An alternative approach is based on the fact that each collision of a photon changes its direction of propagation 
at random. Its motion is thus a random walk. The mean square length of each step is iV^, where i is the mean 
free path (because the distance travelled between two successive collisions obeys Poisson statistics). The theory of 
Brownian motion implies that the mean square distance traveled after n steps is (r^) — 2Pn — 2ict, in agreement 
with the solution of the diffusion equation. The mean free path given by Eq. ( |7^ ) equals £ = 1.8 cm, and the average 
number of collisions of a photon during its travel from the center to the surface is thus n = 8 x 10^". 

C. The temperature within the star 

Consider a stationary regime. The power Q — J (Prq{r) produced in the active core of the star should exactly 
balance the luminosity L. Equations (^-(77) yield 



dr \ pK dr 



(86) 



Taking into account the boundedncss of r(0), Eq. ( p6[ ) provides for a given p{r) and q{r) the temperature profile in 
terms of a single boundary condition, T(Q) or Tg. 

In a model where p is regarded as a constant and where q{r) is also taken as uniform in the core with radius i?c, 
the solution of Eq. (§^) is 
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T\r) = T\0) - {r < R,) (87) 

T\r) = T\0) - ^ (sRi - 2:^) . (r > i?e) (88) 
8cr V r / 



For r — R, Eq. ( p8[ ) relates the surface temperature Tg = T{R) and the central temperature T(0) 



For the Sun, if we take i?c = tsRq and if we make use of Lq — Q — ^nR^q, we find from Eq. (RSh the central 



temperature 



T(0) 



QpnQ 



327rcri?c 



1/4 



SxlO^K. (89) 



This temperature is determined mainly by the heat production rate and by the properties of the core (opacity, density, 
and size). Its insensitivity to the outer regions is related to the efficiency of the photons confinement. The temperature 
of the core is nearly uniform as T{Rc) = 0.9T(0). The value of the temperature in Eq. (|89| ) is not sufficient (by a factor 
of 2) to ensure the ignition of the thermonuclear reactions. The actual value within the Sun, approximately 10^ K, 
can be explained by the fact that we underestimated the opacity k. Actually, besides their Thomson scattering off 
the electrons, photons can interact with partly ionized atoms such as Fe which, although scarce, have large absorption 
cross-sections. Moreover, the non-uniformity of the density enhances the confinement of photons near the center 
where p{r) has a maximum, and hence raises the central temperature above the value of Eq. (p9|) . 



Outside the core and up to the surface, the temperature given by Eq. (88) decreases as 

r^(r) - T^jr) _ 2Rc R-r 
T4(0) ~ 3i? - 2i?c r ■ 



(90) 



D. Relation between mass and luminosity 

The luminosity L of a star is related to its surface temperature through the Stefan-Boltzmann law, Eq. (|l|), and 
to the internal temperature gradient, through the diffusion equation for photons, Eqs. ( [76| ) and ([77|). The temperature 
outside the core thus satisfies 

^_16^^^ (91) 
3pK dr 

Equation (pVt) is consistent with Eq. (p3) and L — Q. 



We have found in Sec. IV the relation ( |38f ) between the mass, the radius, and the internal temperature of a star. 



By eliminating the temperature, we can therefore express the luminosity in terms of the mass and the radius. From 
Eq. ( po| ) we find for —r^dT^/dr the estimate ^RcT'^ifi). The average temperature entering Eq. ( ^8|) is found from 
Eqs. (i^) and to be {R^/ RY'^T{Q). We thus obtain in SI units, using Eq. <M\ and the expression for cr. 



^ ^ _ 32^ f Gm,M V ^ 7.7xlO^^G^m^ ^ l^^^^ 

9pK V lO/si? / h^c'^K K ' ^ ' 



Letting M = Mq and taking for k the value in Eq. (83), we find L ~ 2.2 x 10^® W, in reasonable agreement with the 
actual luminosity Lq = 3.8 x 10^^ W of the Sun. 

Remarkably, the luminosity is found to depend only on the mass of the star, with a power law L cx M'^ which is 
consistent with the observational data for AIq < M < 100 A/q. The fact that the radius (or the temperature) does 
not enter Eq. (^2|) implies that the luminosity L is independent of the power production Q. In the stationary regime 
L and Q are equal, but perturbations around this regime do not modify the luminosity although Q is changed. We 
can recover the result L oc by means of a qualitative argument, using L oc RT'^/ pn cx {RT)'^/Mk, and, from the 
virial equation (|3^), RT cx GM. The coefficient is found within a numerical factor by dimensional analysis. 

For smaller stars where microscopic radiative processes other than Thomson scattering occur, the opacity k is 
not a constant but varies as pT~''/^. It will be left as a problem to show that L cx M^T^/^ for such stars. 

A. Nuclear reactions VII. THE NUCLEAR FURNACE: BALANCE EQUATIONS 
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The source of heating in the core of the stars of the main sequence is thermonuclear fusion. For M < 0.5 Mq the 
dominant process is the sequence of reactions 

p + p ^ d + c+ + 1/ + 0.4 MeV, (93) 

e+ +c- ^ 7+ l.OMeV, (94) 

d + p^ ^He + 7 + 5.5MeV, (95) 

^He + ^He^ *He + 2p+ 12.9 MeV, (96) 



which altogether amount to 



4p + 2e" -> "^He + 2iy + 26.7 MeV. (97) 



The same reactions also take place in the first stage of the life of heavier stars, with 0.5 Mq < M < 100 Mq. However, 
in such stars after these reactions have taken place in the core and after the temperature has sufficiently risen, 
new fusion reactions involving He can occur. Strong interactions are involved in the reactions ( p5| ) and (pq), and 
electromagnetic interactions in the reactions ( p^ and (|95|). These three reactions have large cross-sections compared 
to (|93|), which involves weak interactions, and which is a prerequisite for the other ones to take place. The reactions 



(|93|)V(|95|), and (96) are moreover hindered by the Coulomb repulsion between the two nuclei, which must be overcome 
before they can fuse. 

The reactions (p3|)-(p6|) are all exothermic, and energy is released in the form of kinetic energy of the particles 
produced. Collisions distribute this energy to other particles, and thermalization takes place, so that the effect is 
heat production. However, the neutrinos created by the reaction (|9^) escape the star without having interacted 
with matter, so that their energy, on average of 0.3 MeV per neutrino, is lost. The final balance is thus an average 
production of heat of 77 = 13 MeV per fusion of a pair of protons. 

The cross-section for the reaction (^3|) is especially small because this reaction combines weak interactions with the 
need for the two protons to approach each other at a distance comparable to the range rg of nuclear interactions, of the 
order of 1 fm. For this reason, the reaction ( p3| ) cannot be produced in a laboratory, and its cross-section needs to be 
evaluated theoretically. The value of the Coulomb repulsion at the separation ro — 1 fm is e'^/Aireor = 1.3 MeV; thus, 
even at the temperature of 10^ K (corresponding to kT = 0.9 keV) at the core of the Sun, the classical Boltzmann- 
Gibbs probability of passing over the Coulomb barrier is negligible because exp[— 1.3 x 10^/0.9 x 10'^] ~ 10^^°*^. 
Quantum tunneling through the barrier is thus the only process that allows the reaction (|9^). Its probability is given 
by the WKB approximation, so that the cross-section of a reaction between two charged nuclei has the general form 



S{e) 

exp 



rfZ^Z2e 



1/2 



dr 



(98) 



The relative mass is denoted hy fj. = mim2/(mi + 7712), and the relative energy by e = ^/^(vi — V2)^ in terms of the 
velocities of the colliding nuclei. The factor 1/e has a geometrical origin, and the remaining coefficient S{e) which 
accounts for the short range nuclear interactions varies slowly with energy (except near resonances). The integral 
runs from the range ro of the nuclear forces to the turning point, but the lower bound can be replaced by 0, which 
yields 

= ^exp(-v/i^), £B = ^g^,| . (99) 
The energy £b, which characterizes the strength of the potential barrier, is large (500 keV for the p-p reaction) 



compared to the thermal energy (l keVat 10^ K) . The exponent in Eq. (99) is thus small, and varies rapidly with 



e. For the weak interaction (93), the coefficient S(e) may be estimated by using the Fermi golden rule, the coupling 
constant being inferred from the /3-decay half life of the neutron.^ Its value, S ~ 3.8 x lO"^*^ m^ keV, is smaller than 
for strong interactions by more than 20 orders of magnitude: for the reaction (BO), an extrapolation of data obtained 
by ion collision experiments at higher than thermal energies yields S* ~ 5 x 10^^^ vc? keV. 

B. The Gamow peak 

The probability of a single reaction induced by the collision of two positively charged particles is characterized 
by its cross-section which is of the form (|9^). We wish to evaluate from it the average number w of such reactions 
which take place per unit time in a unit volume where the particles constitute a gas in thermal equilibrium. If particle 
2 with mass and velocity V2 is regarded as a target, the relative flux of particles 1 having the density nx and a 
velocity within the volume d^Vi around vi is, according to Maxwell's distribution. 
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\^^-^^\<^f'"^~"''^^''""^d'^^- (100) 



The corresponding average number of reactions per unit time of the particle 2 is the product of Eqs. ( |99| ) and (100) 
The average number w is then obtained by summing over the particles 2 as 



I mi .3/2. 7712 x3/2 /■ 3 3 , , 



5(e) 
X— ^exp 

£ 



[-^I^e-{mivl + m2vl)/2kTY (101) 



Integration over all variables except e = 5m(vi — V2)^ yields 

4nin2 



JdeS{e)exp{-y^I^-e/kT). (102) 



(27r/i)V2(fcT)3/2 

Due to the rapid increase with e of the tunneling factor in the cross-section, and to the exponential decrease of 



the Boltzmann factor, the integrand in Eq. (102) has a sharp peak, called the Gamov peak. The position eq of this 



peak, refered to as the Gamov energy, is the value of s for which the exponent is largest, namely 

The Gamow energy is much larger than the thermal energy, and much smaller than the barrier energy eb- The width 



of the Gamow peak is of order (fcrec)^^^- The particular shape of the integrand in Eq. (102) implies that most 
reactions take place at a relative energy in the range eq ± (kTeQ)^^^ . Slower particles are not sufficiently energetic 
to pass through the barrier, while the Maxwell distribution cuts off the number of faster particles that would collide 
more efficiently. A compromise is thus needed, which results in the Gamow peak. 

Neglecting the variation of S{e) within the Gamow peak and replacing the exponential in Eq. ( |102| ) by a gaussian, 
we obtain finally 



4(2£ 



|l/6 



" = (3j)V2(fcr)2/3 "^"^ • 

The exponential factor, with sb defined in Eq. (p9|), comes from both the Coulomb barrier and the Maxwell factor, 
because (eB/4fcr)^/3 — (eB/4eG)^^^ = ec/kT. It dominates w, increasing rapidly with the temperature. 

For the first reaction (p3|), we have Zi = Z2 — 1, — ^mp, Eb = 490 keV, and a factor i should be included 



in (104) due to the indistinguishability of the two colliding protons. For the Sun, a numerical estimate of w can be 
obtained by assuming that the density of protons in the core is 100 times larger than the average density (see Fig. 1), 
that is, ni = n2 = 10^^ m"^^ and that the temperature is kT = 1 keV; by using S = 3.8 x 10~^° m^ keV we find a rate 
w = 2x lO^-^m-^s-i. 

C. Povifer production 

Whereas the major part of the energy is generated by the strong interactions (p5| ) and (p6|), the time-rate of 
the entire process is controlled by the interaction (9^) which is both weak and hindered by the Coulomb repulsion. 
As soon as a positron and a deuteron are produced by the reaction (p^, they readily find in their neighborhood an 
electron and a proton to react according to (|^ and (^), which have much larger cross-sections. The final reaction 
( p6| ) also has a large cross-section, although the probability of tunneling is strongly reduced by the fact that eb in 
Eq. ( p9[ ) is 48 times larger for this reaction than for the p-p reaction. However, due to the factor [^(^He)]^ which 
enters Eq. ( |104| ), the reaction (^6|) occurs less often than (^3|) as long as a sufficient amount of ^He has not yet been 
accumulated (which typically takes 10^ years). The stationary regime is attained when n('^IIe) has reached a value 
such that the rate at which ^He is produced by the reactions (|9^) and (^) equals the rate at which it disappears 
through the reaction (pw. 



This stationary regime is therefore governed by the rate w evaluated above for the p-p reactions (93). Because 
each fusion of a pair of protons releases altogether an energy 77 = 13 MeV, the production q{r) of heat per unit time 
and unit volume (or unit mass) is 

q = T]w~ 500Wm"3 ~ 4 X 10"^ Wkg"\ (105) 
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This value of q is small, smaller for example than the heat radiated by a human body, which is of the order of 
1.5 Wkg"""^. This small value implies that the nuclear reactions do not significantly perturb the state of the Sun. The 
total heat power Q is large due to the huge size of the Sun. By assuming that the active core is a sphere of radius 
Rc = i^-R©, we find 

Q = Ijri?^ g ~ 2 X 10^^ W, (106) 

in reasonable agreement with the actual luminosity of the Sun, Lq = 3.8 x 10^^ W. 

The small efficiency of nuclear fusion in stars is illustrated by the average time it takes for a particular proton 
in the core to react, ^nw~^ w 8 x 10^ years. Accordingly, the total time during which the Sun can radiate as it does 
presently is large; assuming that 10% of its mass can be burnt, we find 

V Mq 1 ^^10 



2 lOmp Lq 



10'^ years. (107) 



The flux of energy carried by photons is accompanied by a flux of neutrinos, produced by the reaction (631). According to 
Eq. (p?!), one neutrino should be emitted for each 13 MeV released. Thus the flux of neutrinos emitted by the Sun can be 
evaluated from its luminosity. However a deficit exists in the detection of the solar neutrinos on Earth, which has not yet been 
explained. 

More precise calculations should account for the non-uniformity of the density and the temperature, and for the presence of 
elements other than H. The latter efl'ect is especially important for the heavier stars and for the later stages of evolution, when 
the central part of a star becomes much hotter and much more dense than the outer shells and when elements other than 
H and He have already been synthesized. Indeed, while the onset of the primary cycle (^^-(^6|) of reactions takes place at 
temperatures of the order of lO'^ K, the onset of other reactions requires higher temperatures because the Coulomb barrier in 
Eq. (BQI) increases with the charge and the mass of the colliding nuclei. In red giants with masses Mq < M < 10 M©, H can 
thus fuse into He outside the inert helium core, but through reactions other than the cycle (p3[)-(p6[) which are catalyzed by 
the presence of C. Helium can in turn fuse into C provided the temperature reaches 10* K (instead of 10^ K for H). In this case, 
an additional effect requires a high temperature: whereas the reactions in (p^)-(|9^ are exothermic, the successive reactions 
*He-|- *He-|- 92 keV — > *Be* and *He-l-*Be* -|-67keV — > ^^C* are endothermic, and the thermal energy of ^He should overcome 
this threshold. Heat is thereafter retrieved through the electromagnetic deexcitation ^^C* ^^C + 7 + 7.4 MeV. 



D. The stability of stellar equilibrium 

The stars in the main sequence reach a stationary regime in which their luminosity L equals the heat production 
rate Q: the radiated energy is then exactly compensated for by the thermonuclear power. This stationary state, 
refered to by astrophysicists as "equilibrium," lasts for a long time in the life of the star until the nuclear fuel is 
exhausted. 

The stability of the stationary regime is paradoxical because the production Q of heat is a rapidly increasing 



function of the core temperature, as shown by Eqs. (104) and (|106|), whereas the luminosity L depends practically 



only on the mass, as shown in Sec. VI D. To study perturbations around the stationary state, we use the dynamical 
balance equation: 

dE 

^ = Q-L, (108) 
which expresses energy conservation. If a s mal l perturbation slightly raises the core temperature above its value for 



which Q — L, the right-hand side of Eq. (108) becomes positive, and the total energy E increases. The very slow 
progression of heat from the core to the surface (typically lO''^ years) prevents the extra heat that is produced from 
being radiated. If we naively accept the idea that the total energy E varies in the same direction as the internal 
energy U, the initial perturbation would result in an increase of U, and hence in an extra increase of the temperature. 
The heat production Q would rise again, and so on, leading eventually to a heating instability and to an explosion. 
If conversely the core of the star happens to cool down slightly below the temperature of the stationary regime, the 
thermonuclear reactions would by the same hypothesis stop progressively. 

Fortunately, the total energy E = Eq + U is dominated not by U, but by the gravitational contribution Eq, equal 
to —2U when the star is in gravitational equilibrium. We noted in Sec. pV] C that the time-scale for the establishment of 
gravitational equilibrium is short, of the order of a few hours. Hence, an accidental increase of the core temperature 
Tc which leads to an increases of Q and thus of E, produces an immediate mechanical reaction: the star slightly 
expands according to Eq. (|3^), coohng the core of the star back to its original temperature. Conversely, an accidental 
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lowering of the thermonuclear activity produces a decrease of E, and restoration of gravitational equilibrium results 
in a contraction and a heating so that the rate Q of heat production returns to its stationary value L. The stability 
of the stationary state of a star is therefore due to the negative sign of dE/dT, an unusual feature that we noted in 



Sec. |V|B. 

Thus, the regulation of the thermonuclear power plant constituted by the core of a star is ensured by the 
gravitational forces, which prevent Q from shifting away from L. 

The state of a staVlila gii^BffliHiESiKclfiaeiEBiitl^nMiPqEMQIiliTiBNT^ Sa^ARSnature and proportion 
of the constituent particles; the local mass density p(r) and the radius R; the pressure P{r) and the local energy 
density u{r)] the temperature T(r) and its extreme values at the surface Tg and at the center T^, the total energy E, 
including the gravitation energy Eq and the internal energy U; the local rate q{r) and total rate Q of heat production 
in the core for the active stars of the main sequence; the luminosity L and the spectral constitution of the emitted 
light. The last quantity is the only one that we can observe by photometry of the radiation that we receive. 

In Sees. Ill to VII we have used statistical mechanics to establish relations among all the quantities listed above. 
We will now summarize the foregoing results to show that a single quantity, the mass M of the star, essentially 
determines all of the properties of the object as a function of time. 

Actually, a star always begins as a cloud of hydrogen, with about 25% in mass of helium and with small quantities 
of heavier elements. The presence of these elements may influence some properties, in particular through their 
interactions with photons (Sees. ^ and VI). They may thus play an important role at some stages of the evolution 



of stars, but they do not govern the main common features of stars. On the other hand, we have considered only 
isotropic objects; some peculiar features can emerge for other geometries. Binary stars rotating close to each other 
can in particular produce accretion disks or explosive events (novae). 

For a given constitution of matter, equilibrium statistical mechanics provides the equations of state which relate 
P{r) and u{r) to p{r) and T(r). The ionization ratio and the criterion for using quantum mechanics or relativity 
depend only on p{r) and T{r). 

The constitution of the matter of stars is determined mainly by their age and by their mass which governs 
their evolution. Indeed, like protostars, stars of the main sequence are mainly made of hydrogen plus an increasing 
amount of helium and of other elements produced by nuclear fusion. White dwarfs are made of the latter elements, 
in proportions depending on the degree to which fusion has been effective, itself determined by the temperature of 
the core which is hotter for heavier stars. For neutron stars produced in the supernovae explosions, nuclear fusion is 
complete. 

The gravitational energy Eq results from the density through Eq. (|l^) and is related to the pressure through the 
virial equation (|50|). The hydrostatic equilibrium equation ( ^8|) provides the function p{r), and hence the radius R 
and the total energy E, in terms of the function T{r) and the mass M. 

The temperature T{r) is in turn a solution of the heat transport equation (p6|), which relates the surface temper- 
ature Tg to that of the center Tc and to the rate q{r) of heat production. This solution depends on a single parameter, 
say Tc. 

The luminosity L of the star then follows from the surface temperature. The more detailed characteristics of its 
radiation can be obtained by solving the Boltzmann-Lorentz equation (^) in the outer shells. 

The theory of thermonuclear reactions expresses the heat production Q as a function of the temperature and the 



density in the active core through Eq. (104) 



All the above quantities are thus related by means of a set of self-consistent equations to the mass M and to the 



only remaining parameter Tc. Finally, conservation of energy yields the differential equation (108), which determines 
the time-dependence of the parameter Tc and hence of all quantities. Altogether, apart from minor differences, all 
stars having the same mass have the same history. 

The classification and evolution of stars according to their mass can be visualized by the diagrams of Figs. 1 
and 2 or the schematic diagram of Fig. 3. As an exercise, the various features of Fig. 3 can be derived from the 



results of Sees. Ill and IV. For protostars, white dwarfs and neutron stars, we have Q = in Eq. ( |10§| ) and energy is 
gradually lost. The evolution is represented by constant mass lines in the various figures, except for supernovae. The 
time-dependence is not uniform. In particular the evolution stops for a long time for the stars of the main sequence. 
As a general rule, more massive stars evolve faster. This behavior can be illustrated by the duration of the stationary 
regime of a star in the main sequence, the luminosity of which is proportional to M^. Because the quantity of nuclear 
fuel which supplies this radiation is proportional to M, the lifetime of the star is proportional to 1/AP. Massive stars 
thus waste their nuclear stocks. 

The simple model that we often used, neglecting the variations with r of most quantities, is adequate in many situations, but 
some stars need to be analyzed in shells having different qualitative properties. For instance, the study of heat transfers in 
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white dwarfs (Ref. 2, Exer. 15f) requires that we distinguish the interior, which is a completely ionized gas of elements such 
as C or O, from the envelope, which is a colder gas of non-ionized, lighter atoms. Heat conduction takes place in the interior 
through transport and collisions of electrons on the nuclei like in a plasma or a metal, and it is so efficient that the interior 
temperature is uniform. In the envelope the transport takes place through radiative processes as in Sec. IVTO, which shows that 
the surface temperature is anproximately proportional to the internal temperature, Ts ~ r/2000. Noting that the specific heat 
includes the contribution (|l7| ) of the electrons and the classical contribution of the nuclei, we can then study the duration of 
the glare of white dwarfs. Starting from an initial temperature T of 2 x 10^ K, it is found that T decreases very slowly, by 10% 
in 2 X 10^ years. The statistics of the observed luminosities of such stars thus reflects the statistics of their ages. It is consistent 
with a constant rate of creation of white dwarfs in our galaxy — except for the occurrence of a clear cut-off toward the faintest 
objects: no white dwarfs are seen with luminosities below 2.5 x 10~^Lq. This anomalous absence of very old, cooled white 
dwarfs can be used to estimate the age of our galaxy, which is thus found to be of the order of lO^'^ years. 

Another example of heterogeneous stars are the red giants, a stage reached by the stars of the main sequence with masses 
Mq < M < 10 Mq after hydrogen in the central part has been converted into helium. These stars have a large radius, typically 
100 Rq , and hence a large luminosity although their surface temperature is lower than that of the Sun. These properties can be 
understood by the use of a two-zone model, including a very dense and hot inert core of helium, and a huge hydrogen envelope, 
the deepest layer of which is heated by the gravitational contraction of the core and undergoes thermonuclear reactions. 

For very bright stars, the contribution of photons to the energy and pressure becomes important. The stability of the dilute 
envelope of such a star can be studied as a problem. According to Eqs. ( p^ ) and ([77[), the radiation pressure P-y satisfies 

dr AiTcr'' 

where Konv is the opacity of the envelope. Typically Kenv = 5ft;, where k is the opacity in the bulk. From Eq. (^sj), we infer 
the existence of a critical luminosity AncGM / n^m above which the envelope is expelled by radiation. An upper bound for the 
mass of stars is thus obtained from the mass- luminosity relation (|92|), 



M„,ax = ^ T7 . (110) 



G \k, 



which yields Mmax = 100 Mq . 

We have considered stationary or slow regimes only, in which the evolution is due to radiation and governed by Eq. (108). 
More rapid dynamical processes can suggest additional student problems. The above equations should then be supplemented 
with conservation equations, involving the time-derivatives of the energy density it and the photon density for transport of 
heat or of radiative energy as in Sec. VI, and time-derivatives of the mass and momentum densities for transport of matter. 
These time-derivatives are related to the corresponding thermodynamic forces. When local equilibrium is not reached, one 
should instead resort to kinetic equations as in Sec. ^ Various phenomena can be understood by solving the resulting partial 
differential equations. For instance, simple models explain the solar oscillations. Thermal convection, which also exists in the 
Sun, governs the transfers of matter and heat occurring in the envelope of red giants. Important quantities of matter can be 
expelled at some stages of evolution, in particular during the transformation of a red giant into a white dwarf, and especially 
in heavy stars, which for M > 40 Mq can loose half of their mass due to radiation pressure. Dynamics is essential in all 
the transient regimes. For instance, when the ignition temperature of helium is reached in the core of a red giant, a fiash 
occurs, which rapidly brings the star from one stationary regime to another. When nuclear fuel is exhausted, processes are also 
rapid because nothing opposes the gravitational contraction. Depending on the mass, the resulting instability may produce 
pulsations, implosions, or explosions. The dynamics and luminosity of accretion disks and the understanding of novae provide 
problems of thermodynamics and statistical mechanics accessible to students. 

We have successfully proposed to students a series of exercises explaining the main features of supernovae and involving only the 
above ideas. We now sketch what can be learned from these exercises. In its preliminary stage, a supernova is a red supergiant 
with typical mass 10 Mq and radius 5 x 10* km, made of a compact core of ^^Fe at a temperature of 5 x 10^ K surrounded 
by successive shells made of ligh ter and lighter elements. The mass of the core grows slowly from 1 Mq by nuclear fusion. It 



is shown as in Sees. Ill and |V| that the thermodynamics and the gravitational equilibrium of the core are dominated by the 
ultrarelativistic Fermi gas of electrons. Hence, beyond some mass analogous to the Chandrasekhar limit mass for white dwarfs, 
the core no longer tends to expand, being maintained by the pressure of the envelope, but tends to shrink. This gravitational 
instability produces a sudden collapse when the critical mass is reached. The duration of the collapse can be estimated as the 
time of free fall of a particle under the effect of gravity from a height equal to the initial radius of the core, say 2000 km, which 
yields a fraction of a second. The collapse produces heating and fusion of the Fe nuclei into a unique gigantic single nucleus, 
which is a neutron star. Stability is then recovered because the equation of state becomes that of the non-relativistic Fermi gas 
of neutrons. Most of the energy of this implosion is released through neutrinos formed in the nuclear reaction, which can escape 
freely. Their number can be estimated by assuming their average energy to be 10 MeV and by comparing the energies of the 
initial Fe core and of the resulting neutron star (a difference of 10^^ J is found). One can also explain the number of neutrinos 
detected on Earth during the collapse of the supernova 1987A, which is 150 000 light-years from us. Part of the energy produces 
a shock wave in the envelope which explodes, in the same way as the shock on the ground of a pair of elastic balls falling 
together on top of each other expels the upper one to a high altitude (an experiment easily done in class). Assuming that, 
after the shock wave has passed across the envelope, the kinetic and thermal energies of matter are equal, one can estimate the 
speed of expansion and the temperature of the envelope. As a consequence, one can show that the energy is mainly carried by 
the photons in thermal equilibrium, and hence that the temperature varies as 1/T during the subsequent rapid expansion. The 
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luminosity after an initial intense flash can then be evaluated as in Sec. IVJ. ft is found to remain nearly constant in time for a 
few months during which an isentropic expansion takes place. Its value, oForder 10^ L©, and the duration of the flare-up can be 
estimated by assuming that the shock-wave carries 1% of the total energy released, that is, 10"*^ J. Afterwards, when the total 
energy radiated becomes equal to the residual thermal energy, the expansion is no longer isentropic and the luminosity rapidly 
decreases while the cloud formed by the envelope goes on expanding. It is found that the total energy radiated by photons is 
1% of the kinetic energy, and 10~* of the energy is carried away by the neutrinos. 
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Figure Captions 

Figure 1. Protostars and stars in the density-temperature plane, where the coordinates are logrie and logT in 
units of and K, respectively. The curves describe the evolution for a fixed mass M (in units of M©), obtained 
from Eqs. ( ^o|) and ( ]4l| ) with = A^nuc- While the density increases, the temperature rises, proportionally to n^^^ in 
the classical regime, and then decreases. The maximum temperature (dashed line) is given by Eq. (p4[). The dotted 
line indicates the ionization boundary of Eq. ( |l5| ) in the classical regime. When the mass is sufficient (M > 0.1 Mq), 
nuclear ignition takes place and the evolution of n and T stops. The star indicates the stationary state of the Sun. 
Note that its coordinates represent average density and temperature; the actual temperature of the core is higher (10^ 
K). White dwarfs are obtained after exhaustion of their nuclear fuel and nearly complete contraction in the extreme 
quantum regime. 

Figure 2. Protostars and stars in the density-energy plane, where the coordinates are logrio and log|_E|/iVo = 
log |P/nc in units of m^'^ and eV, respectively. The same conventions as in Fig. 1 are used. The thin curves are the 
isotherms (labelled by logT in K): in the classical regime they are constants, reflecting the equipartition theorem, and 
they bend up in the quantum regime. The entire evolution, contraction, and loss of energy of a star with mass M (in 
units of Mq) is represented by a straight line with slope 1/3 ending up on the dashed-dotted line which marks the end 
of the evolution (Jupiter- like objects and white dwarfs). The dashed line is the locus of the maximum temperatures 
reached during the evolution. 

Figure 3. Evolution of the stars in the mass-density plane with coordinates log M /Mq and log p [p in kg m~^) . The 
evolution takes place upward, with a loss of mass for massive stars when matter is expelled after exhaustion of hydrogen 
and overheating. The thin lines represent the isotherms of Eq. (pq), for average temperatures corresponding to the 
ignition of H (core temperature 10^ K), He (10* K), C and O (10^ K); we assume the core temperature to be 10 times 
the average temperature. The stars in the main sequence correspond to the isotherm 10^ K for 0.1 Mq < M < 100 Mq. 
The isotherms bend up for large masses as an effect of the radiation pressure; for small masses they bend down (not 
drawn) toward the line of dwarfs due to the Fermi gas behavior of electrons. The dotted line is the ionization boundary 
of Eq. ( [Tsl ) . The dashed-dotted lines mark the end of the evolution for aborted stars and white dwarfs (with M less 



than the Chandrasekhar mass 1.4 M©), and for neutron stars (Eq. (P3|)). The maximum temperature (44) is attained 
along the dashed line. The thick dashed line is the black-hole boundary of Eq. (E6h 
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